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ABSTRACT 


Specular  returns  from  radar  targets  contain  large  power  which  makes  a 
technique  based  on  the  analysis  of  specular s  attractive.  Herein,  the  specular 
returns  from  a  flat  region  and  from  a  smooth  surface  tangent  to  a  plane  along 
a  curve,  which  are  the  largest  returns  for  reasonably  sized  bodies,  are  used 
to  obtain  information  about  the  geometry  of  the  target. 
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INTRODUCTION 


Specular  returns  from  scattering  targets  may  be  classified  according  to 
the  order  of  magnitude  of  the  return.  It  can  be  shown  that  some  typical  returns 
are  (after  suitable  normalization)  of  the  following  order  in  ka,  where  k  is 
the  wave-number  and  a  is  a  typical  target  dimension: 

Order  Surface  Geometry 

flat  region 

smooth  surface  tangent  to  a 
plane  along  a  curve 

Convex  surface  with  non-zero 
Gaussian  curvature 

edge-diffraction 

vertex-diffraction 

Because  of  the  large  power  of  a  specular  return  it  would  be  convenient  to 
use  it  to  obtain  information  about  the  geometry  of  the  corresponding  specular 
region.  In  this  paper  we  develop  a  method  for  doing  this  for  the  first  two  types 
of  surface  regions  in  the  above  table. 


0  [(ka)2] 

o  U3/2] 

0  [(ka)] 

0  [(ka,1/2] 
0  [l] 


QUASI-MONOSTATIC  SCATTERING 


Let  be  a  unit  vector  in  the  direction  of  propagation  of  the  plane  wave 


ikl  -  X 
u.  =  e - 


(1) 
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which  is  incident  on  a  target.  If  the  total  (scalar)  field  satisfies  the  boundary 
condition  u  =  o  on  the  surface  of  the  target,  then  according  to  the  Kirchhoff 
(physical  optics)  approximation,  the  scattered  far-field  is  given  by  (see  Equa¬ 
tion  (2.  25)  of  [1]  ) 


s 


.  ikr 

1  e 

2  :rkr 


(2) 


where 


g 


-  -kH  ±  * 
£ 


ik  (I  -  £)  •  X 


ds(X). 


(3) 


Here  k  is  the  wave-number,  r  is  the  range.  N  is  the  outward  unit  normal 
vector  on  the  surface  of  the  target,  J  is  a  unit  vector  pointing  from  the  target 
to  the  point  of  observation,  and  the  surface  integral  of  Equation  (3)  is  taken 
OA'er  the  illuminated  portion  £  of  the  target.  A  similar  formula  holds  for  the 
boundary  condition  9  u/9  n  =  o  ,  and  for  scattering  of  an  electromagnetic  wave 
by  a  perfectly  conducting  target  (see  [l]  ). 

A  rectangular  coordinate  system  such  that  I_  =  (0,  0.  -1),  J  = 

(sin  9,0,  cos  9),  and  X  =  (x,  y,  z)  is  introduced.  In  the  quasi-monostatic 
case  9  «1  and  J  ~  (9  .  0  ,  1).  Then  since  -T  '  N  ds  =  dxdy.  Equation  (3) 
becomes 

....  i  2  Cf  -ik[9x  +  2z(x,  y)]  ,  ,  . .. 

g  (9)  =  k  J  j  e  ^  v  yn  dxdy  ,  (4) 

P 

where  P  is  the  projection  of  £  onto  the  xy-plane. 

[l]  R.  M.  Lewis,  Notes  on  the  Kirchhoff  Method  in  Scattering  Problems. 
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For  simplicity  it  is  now  assumed  that  the  region  P  is  intersected,  at 
most,  twice  by  each  line  parallel  to  the  y-axis.  Then  its  boundary  is  given  by 
the  two  functions  y  (x)  and  y9  (x)  as  illustrated  in  Figure  1.  These  functions 
are  defined  so  that  they  are  continuous  and  have  constant  values  outside  of  P. 
(Then  y  =  y  outside  of  P).  We  take  2L  greater  than  the  largest  diameter 

J.  Z 

of  P  parallel  to  the  x-axis  and  choose  the  origin  so  that  the  points  ±  L  lie 
outside  of  P. 

From  Equation  (4)  we  now  see  that 

L 

C  .  .  -lkQx 

J  y  (x)e  dx,  (5) 

-L 

y2W 

£  e-2U“(x'y)dy.  (6) 

yl 

From  the  theory  of  Fourier  series  it  is  known  that  if  f(x)  is  piecewise  smooth 
and  if  the  Fourier  coefficients 


g(9)  =  k 


where 


y(x)  =  k 


a 

n 


L 


-L 


f(x)  e 


-in7rx 

L 


(7) 
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are  introduced,  then 

oo 

f(x)  =  ^  anelT 

_  n  L 

OO 

n  =  -oo 


Therefore,  if 


is  set,  then 


where 


mr_ 

kL 


o,  ±1  ,  ±2, 


Y  (x) 


C  e 
n 


iriTrx 

L 


n  =  -°o 


C 

n 


— 

2kL  s  kL 


(8) 


(9) 


(10) 


(11) 
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Thus  if  g  can  be  measured  for  those  values  of  9  given  by  Equation  (9) 
y  (x)  can  be  determined.  The  latter  function  is  closely  related  to  the  geometry 
of  the  target  as  illustrated  in  the  next  section.  A  multi- static  system  for 
making  the  required  measurements  instantaneously  is  illustrated  in  Figure  2 
(other  measuring  systems  can  obviously  be  used). 

As  illustrated  in  Figure  2,  receivers  are  placed  on  a  straight  line  at 
distances  from  the  transmitter  that  are  an  integral  multiple  of  t  .  The  system 
is  characterized  by  four  length  parameters: 


r  range 

t  receiver  spacing 

2L  body  dimension  (upper  bound) 

X  =  2ir/k  wave-length 

Since  r  9  ~  n  t  ,  we  see  from  Equation  (9)  that  the  four  parameters  must 
satisfy  the  condition 


r  _  X 
r  2L 

In  practice  one  would  first  choose  2L  larger  than  the  maximum  expected 
dimension  of  the  target  and  then  set  X  =  2Lr/r. 


(12) 


Of  course  it  is  impractical  to  measure  except  for  small  n,  say 
n  =  o  ,  ±1,  ±2.  Therefore,  it  is  envisioned  that  measurements  will  be  made 
only  of  specular  returns.  For  such  returns  the  larger  Fourier  coefficients 
will  be  negligible  and  a  small  number  of  terms  in  Equation  (10)  should  give  a 
very  good  approximation  to  y  (x).  This  is  certainly  true  if  y  (x)  is  smooth. 

If  not,  then  a  few  terms  will  yield  a  smoothed  approximation  to  y  (x). 
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Figure  2.  Illustration  of  a  Multistatic  System  for  Making 
Required  Measurements  Instantaneously 


SPECULAR  ANALYSIS  OF  A  FLAT  REGION 


The  largest  specular  returns  are  produced  by  flat  portions  of  a  target. 
Assume  that  the  region  P  in  Figure  1  is  a  flat  surface  region  normal  to  the 
direction  of  incidence  L  Then  one  may  take  z(x,y)  =  0  in  Equation  (6).  Now 


y  (x)  =  kw  (x) 


where  w  is  the  width  function 


w  (x)  =  y2  (x)  -y1  (x)  . 


(13) 

(14) 


Furthermore 


g  (0)  =  k  ^  y  (x)  dx  =  k“A, 

-L 

where  A  is  the  area  of  P.  Thus  the  monostatic  specular  return  of  a  flat 


(15) 
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region  is  of  order  (ka)  as  noted  in  the  introduction.  Furthermore,  since 

y  (x)  is  real,  one  sees  from  Equation  (5)  that  g*(0)  =  g(-9)  .  Hence 

C  =  C  *  .  Thus,  for  a  flat  plate,  the  number  of  receiver  sites  can  be 
-n  n 

reduced  by  a  factor  of  two. 

The  function  w(x)  does  not  uniquely  determine  the  shape  of  the  region 
P,  i.  e.  ,  the  functions  y1(x)  and  y  (x),  but  we  will  see  in  the  next  section 

that  if  a  second  measurement  is  made  corresponding  to  a  rotation  in  the 
xy -plane,  the  shape  of  the  region  P  can  be  found. 


SOLUTION  OF  THE  WIDTH  PROBLEM 


Let  P  be  a  region  in  the  xy-plane  whose  boundary  is  intersected,  at 
most,  twice  by  each  line  parallel  to  the  y-axis,  as  illustrated  in  Figure  1. 

We  introduce  a  rotation 

|  =  x  cos  a  +  y  sin  a,  rj  =  -x  sin  a  +  y  cos  a. .  (16) 

Then 


y  =  |  sin  a  +  77  cos  a  ,  x  =  £  cos  a-77  sin  a.  (17) 

The  rotated  axes  are  illustrated  in  Figure  3. 


Let  x.  (£  ,  a)  and  x  (£  ,  o')  be  the  x-  coordinates  of  the  two  boundary 

X 


points  on  the  dashed  line  |  =  constant  parallel  to  the  r\  -  axis  in  Figure  3. 
From  Equation  (16)  we  see  that  the  functions  x^,a  )are  defined  implicitly  by 


£  =  x.  cos  a  +  y.  (x.) 
J  ‘  J  J 


sin  o; ,  j  =  1,  2 


(18) 


Let  w  (|  ,  of)  be  the  length  of  the  dashed  line  joining  the  two  boundary  points, 
i.  e. ,  the  width  of  the  region  P  in  the  direction  of  the  ri  -  axis.  Then  from 
Equation  (16) 
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V  .  (  X  ) 


It 

V0 

CsJ 

I 

< 


Figure  3.  Rotated  Axes 


W  (£,  a)  =  v2  -  V1  =  [-  x2  sin  o:  +  y2  (x£)  cos  a]  -  [-x1  sin  a  +  yx  (x^  cos  a] 


=  (x1  -  x  )  sin  a  +  [y  ( x  )  -  (x^  ]  cos  a  . 


(19) 


The  following  width  problem  is  considered:  Given  w  (£  ,  a  )  for  a  =  0 


and  a  second  value  of  a  to  determine  the  functions  of  y  (x)  and  y2  (x)  . 


Actually  the  second  value  of  o:  very  near  a  =  o  will  be  taken.  More 
precisely,  the  function  o>^(£,  0)  will  be  used. 

At  a  =  0  ,  we  see  from  Equation  (18)  that  x^  =  £  ,  and  by  differenti¬ 
ating  Equation  (18)  with  respect  to  a  we  obtain 


9  x. 
1 


(20) 
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It  follows  from  Equation  (19)  that  at  a  -  0 

w(x,  0)  =  y2  (x)  -  y1  (x)  (21) 

and  from  Equations  (4)  and  (5)  one  finds  that 

wff(^  .  °)  =  yj  (0  [-y2  (€)]  -y{  (S)  [-y1  ($)]  =  ~  \  “  yi  J  •  (22) 

Since  (see  Figure  3) 


y2(Xo>  =  yl  (xo> 


(23) 


Equation  (7)  may  be  integrated  to  obtain 


(y2  +  yx)  (y2  -  y^  =  y2  (x) 


yx  (x)  =  q(x) 


where 


q(x) 


=  -2  \  w  (£  ,  o)  d£  . 
J  a 


(24) 


(25) 


From  Equations  (21)  and  (24)  one  sees  that 

y2  (x)  +  yj  (x)  =  . 


(26) 


and  Equations  (21)  and  (26)  are  easily  solved  to  yield 


y,  /o  (x)  =  ~  [  -  +  w  (x,  o)]  . 

Jl/2  '  ’  2  L  w(x,  0)  '  n 

Thus  Equations  (25)  and  (27)  provide  a  solution  to  the  width  problem. 


(27) 
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SPECULAR  ANALYSIS  OF  A  SMOOTH  SURFACE  WHICH  IS  TANGENT  TO  A 
PLANE  ALONG  A  CURVE 


As  a  second  example  a  regular  surface  which  is  tangent  to  the  plane 

* 

z  =  o  along  a  curve  y  =  yQ(x)  is  considered.  Important  examples  are 
cylinders  and  cones;  in  these  cases  the  curve  is  a  straight  line.  The  integral 
Equation  (6)  is  evaluated  by  the  method  of  stationary  phase.  A  stationary  point 
occurs  at  y  =  yQ(x)  .  If  the  phase  function  <f>  (y)  =  -2  z(x,  y)  is  introduced, 
then  at  the  stationary  point  $  =  0  ,  <£'  =  0  ,  and  <f>"  =  2z  (x,  y  )  >  0. 


Hence 


.  ,  .  ..1/2  iir/4r 

y  (x)  ~  (7rk)  e  [-z  (x, 


yo)] 


-1/2 


(28) 


But  -z  is  the  curvature  of  the  normal  section  of  the  surface  in  the  direction 

yy 

of  the  y-axis.  If  we  introduce  the  principal  curvatures  K  =  0  and 

2 

K  =  K  =  1/p  ,  then  it  can  be  shown  that  -z  =  K  cos  B  where 

2  yy 

tan  B  ~  y  '  (x)  (see  Figure  4). 
o 

It  follows  that 


y  (X)  ~  e17r//4  V7r kp[x, y  (x)]  [1  +  (y')2]  . 


(29) 


=  J~l 


Since  the  element  of  arc-length  is  ds 

L 

y  (x)  dx  ~  ei7f  k  n/tt  ^  \!~kp 

-L 


+  (y* )  dx  it  follows  that 
XJo' 


g(0)  =  k 


I 


ds  . 


(30) 


*  It  can  be  shown  (by  using  the  formula  of  Rodrigues)  that  the  Gaussian  curva¬ 
ture  vanishes  at  each  point  of  the  curve  and  that  the  curve  is  a  principal  curve 
corresponding  to  the  principal  curvature  K  =  o.  Conversely  if  K  =  o 
along  a  principal  curve,  then  that  curve  is  a  plane  curve. 
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Figure  4.  Smooth  Surface  Tangent  to  a  Plane  Along  a  Curve 


3/2 

Thus  the  monostotic  specular  return  is  of  order  (ka)  as  noted  in  the  intro¬ 
duction. 

We  see  from  Equation  (29)  that  y  (x)  depends  on  the  curve  of  tangency 
y  =  yQ(x)  and  the  non-zero  principal  radius  of  curvature  p  along  this  curve. 

These  two  functions  characterize  the  geometry  of  the  specular  region  in  this 
case  just  as  the  two  functions  y  (x)  and  y  (x)  determine  the  shape  of  the 

flat  region  on  pages  6  &  7.  On  pages  7-9  it  was  proved  how  to  determine  y.^ 

and  y  ,  given  y  .  The  analogous  problem  for  the  determination  of  y  and 

p  ,  given  y  ,  remains  to  be  considered. 
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